Andreev reflection spectroscopy on the heavy- fermion superconductor PrOs4Sbi2: 
Microscopic evidence for multiple superconducting order parameters 
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Point-contact Andreev reflection spectroscopy was performed on single crystals of the heavy- 
fermion superconductor PrOs4Sbi2, down to 90mK and up to STesla. The conductance spectra 
showed multiple structures, including zero-bias peaks and spectral dips, which are interpreted as 
signatures of multiple pairing gaps, one clearly having nodes. Samples with 2% Ru replacing Os were 
also measured, showing pronounced spectral humps consistent with the emergence of a non-nodal 
gap. Detailed analysis of the spectral evolution revealed a magnetic fleld- us. -temperature phase 
diagram for PrOs4Sbi2 characterized by two distinct superconducting order parameters. 

PACS numbers: 74.70.Tx, 74.45.-l-c, 74.25.Dw, 74.20.Rp 



Superconductivity in heavy-fermion materials has been 
a heavily researched topic, particularly for the role that 
strongly-correlated electrons play in the pairing process. 
The discovery of superconductivity in the filled skutteru- 
dite PrOs4Sbi2 has attracted much attention, because 
it is the first heavy-fermion superconductor containing 
neither Ce nor U but Pr atoms which show no ground- 
state magnetic order [l[ . Various unconventional proper- 
ties have been reported in the superconducting state of 
PrOs4Sbi2 [il- Of particular interest is the appearance of 
two superconducting transition temperatures (Tc's) 0,0] 
and low-energy quasiparticle excitations j both phe- 
nomena suggesting an exotic pairing mechanism. Fur- 
thermore, angular magneto-thermal conductivity mea- 
surements have indicated that the pairing symmetry of 
PrOs4Sbi2 undergoes a phase transition in a magnetic 
field 8]. Such a complex order-parameter (OP) phase 
diagram suggests the presence of multiple superconduct- 
ing OPs, reminiscent of the case of UPta @. 

There have been conflicting experimental reports on 
the superconducting gap topology of PrOs4Sbi2, some in- 
dicating the presence of gap nodes while others indicating 
the Fermi surface to be fully gapped [l^]- To reconcile 
the difference between these reports, a recent proposal 
has invoked the multi-band dispersion of PrOs4Sbi2 11 1 
to suggest that there are two superconducting gaps, i.e. a 
nodal gap in a light-mass band and a non-nodal gap in a 
heavy-mass band 12, 13]. Double superconducting gaps 
have been observed in two separate thermal-conductivity 
measurements 14, [l5| . However, whereas the earlier 
measurement indicated both gaps to have s-wave sym- 
metry [13], the more recent measurement indicated only 
one s-wave gap with the other gap having nodal sym- 
metry The latter result suggests that there are 
two distinct OPs with different pairing symmetries and 
Tc's in PrOs4Sbi2 [16]. Although still under debate, this 
multi-band and multi-symmetry scenario could also ex- 
plain the two Tc's seen in heat-capacity measurements on 
Pr(0si_2:Rux)4Sbi2, where the nodal gap becomes dom- 



inated by a no n- nodal gap for Ru-doping above ?sl% 0. 

Despite these prior studies, it is still not clear whether 
PrOs4Sbi2 has one or more superconducting gaps, and 
whether these gaps have the same symmetry. The vast 
majority of these studies have involved hulk probes which 
could not directly detect the effects of sample inhomo- 
geneity. In this paper, we report on point-contact An- 
dreev reflection spectroscopy (PCARS) measurements of 
PrOs4Sbi2 single crystals. PCARS is an inherently lo- 
cal probe and also sensitive to phase of the supercon- 
ducting OP [3]. It has been used to probe the pair- 
ing state of several unconventional su perc onductors , in- 
cluding cuprates [l^, ruthenates 2^ 21], borocarbides 
[2^ , MgB2 [il] and heavy-fermion metals [H, HI]. Our 
measurements on PrOs4Sbi2 showed multiple spectral 
features which provide microscopic evidence for two su- 
perconducting gaps, one clearly having nodes. We also 
measured crystals with 2% Ru-doping, which showed the 
emergence of a non-nodal superconducting gap. Detailed 
analysis of the spectral evolution down to 90mK and up 
to STesla yielded a magnetic field- ws. -temperature phase 
diagram for PrOs4Sbi2, characterized by two distinct su- 
perconducting OPs with different symmetries. 

The single crystals used in our experiment were grown 
with a molten metal-flux method [26j. For the undoped 
PrOs4Sbi2 crystals, electrical resistivity p measurements 
showed a single T(. at w 1.85K with k, 5mK transition 
width, and an upper critical field Hc2 ~ 2.25T at 200mK. 
AC susceptibility x{T) showed two T^'s, i.e. T^i ~ 1.85K 
and Tc2 ~ 1.65K, similar to previously reported results 
[^jli^. PCARS measurements were performed in a '^He- 
*He dilution refrigerator, using Pt-Ir tips on c-axis faces 
of the crystals. The crystals were etched in a 1:1 HNO3- 
HCl mixture to remove any residual Sb flux. Junction 
impedances were typically 0.5-1 $7, resulting in point con- 
tacts which were in the ballistic regime (28| . To minimize 
Joule heating, a pulsed-signal four-point technique [3] 
was used to acquire current- vs. -voltage I-V curves, which 
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were numerically differentiated to yield dl/dV spectra. 
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FIG. 1: (color online), (a) Differential conductance as a func- 
tion temperature for a Pt-Ir/PrOs4Sbi2 point-contact junc- 
tion in zero magnetic field, (b) Temperature dependence of 
various spectral features, after normalizing relative to their 
maximum values. The triangles correspond to zero-bias peak 
height, the squares and diamonds to Si and 52. The circles 
correspond to the excess spectral defined in the text. 

Figure [IJa) shows the temperature evolution of dl/dV 
spectra taken on PrOs4Sbi2 in zero magnetic field. At 
90mK, well below Td « 1.85K, there is a pronounced 
zero-bias peak (ZBP) accompanied by a symmetric dip 
which begins at «±0.4mV, as indicated by 6i in the fig- 
ure. An additional dip is also present at «±0.2mV, as 
indicated by 62- With increasing temperature, the ZBP 
lowers in height and the dips move inward. The spectrum 
flattens out above Td « 1.85K, consistent with the sam- 
ple no longer being superconducting. Figure [TJb) shows 
detailed temperature evolution of the ZBP height and of 
both Si and 62 positions, after normalizing them relative 
to their maximum values. Also plotted in FigUJb) is the 
excess spectral area, which is defined by subtracting each 
dl /dV spectrum by the normal-state spectrum and then 
numerically integrating between ± 1.8mV. These plots 
indicate that the evolution of the various spectral fea- 
tures is governed by two distinct temperature scales, cor- 
responding respectively to Td ~ 1.85K and Tc2 ~ 1.65K. 



The ZBPs seen in our dl/dV spectra can be interpreted 
as distinctive signature of a nodal gap, while the excess 
spectral area provides strong evidence for a non-nodal 
second gap. First, several bulk experiments have already 
indicated PrOs4Sbi2 to have point nodes in its supercon- 
ducting gap structure 0, H, d, [13] ■ Second, theoretical 
modeling of Andreev refiection for PrOs4Sbi2 has shown 
that points nodes in a triplet pair potential can produce 
ZBPs 31|jSimilar to the effect of d-wave line nodes 
for cuprates [34|. Here we point out that ZBPs are a 
manifestation of nodal Andreev states, which are due to 
interference between quasiparticles and thus conserved, 
in contrast to the non-nodal Andreev states, which are 
due to retroflection from pairs and thus not conserved 
Therefore, the excess spectral area seen in FigH] 



indicates that, in addition to contributions from Andreev 
states due to a nodal gap, there are also Andreev states 
due to a non-nodal gap j34| . In fact, since the ZBP van- 
ishes near Tc2 ~ 1.65K while the excess area persists up 
to Tci « 1.85K, these spectral features can be attributed 
to two distinct superconducting OPs In essence, our 
data provides microscopic corroboration for the thermal 
conductivity results of Ref. ^Sj , indicating two supercon- 
ducting gaps with different symmetries in PrOs4Sbi2. 

Figure [H^a) shows the magnetic-field evolution of 
dl/dV spectra measured at 90mK. As the field is in- 
creased, the ZBP lowers in height and the dip structures 
move inward, qualitatively similar to the spectral evo- 
lution with temperature. However, the ZBP height col- 
lapses more quickly with field than with temperature, 
relative to the evolution of the excess spectral area. As 
shown in Figure (Hb), the ZBP height vanishes above 
«1.5T while the excess area persists up to «2.3T. Using 
the two-OP Andreev reflection scenario outlined above, 
this field evolution indicates that there are also two differ- 
ent energy scales governing the field suppression of super- 
conductivity in PrOs4Sbi2. Namely, the nodal Andreev 
states are spectrally robust up to a lower-field bound- 
ary H' « I.5T at 90mK, while the non-nodal Andreev 
states persist up to a higher-field boundary H" « 2.3T 
at 90mK. This two-boundary field evolution is qualita- 
tively similar to the angular magneto-thermal conduc- 
tivity data of Ref. [3], which indicated a phase-diagram 
boundary H*{T) across which the pairing symmetry 
changes. It is worth noting that PCARS data taken on 
UPta which is known to have multiple supercon- 

ducting OPs [§], also shows spectral features that vanish 
at a field boundary lower than i?c2- 

To map out a detailed H-T phase diagram for 
PrOs4Sbi2, we carried out a combined analysis of the 
temperature and field evolutions of our data. Figure [3] 
shows a master plot of the full temperature dependences 
of H' and H" determined from our spectroscopy data, 
along with Hc2{T) determined from our p{T,H) data. 
First, it is clear that H"{T) coincides with Hc2{T), in- 
dicating that the higher-field boundary is just the re- 
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FIG. 3: (color online). Magnetic field vs. temperature phase 
diagram determined from our point-contact spectroscopy and 
resistivity data on PrOs4Sbi2. H' (triangles) correspond to 
the field at which the ZBP vanishes. H" (circles) correspond 
to the field at which the excess spectral area vanishes. Hc2 
(diamonds) is the upper-critical field measured by resistiv- 
ity. Dotted curves are added to guide the eye. It should 
be emphasized that this phase diagram, which indicates two 
distinct superconducting OPs with different symmetries, is 
determined with a highly local probe, thus arguing strongly 
against sample inhomogeneity as the cause of multiple Tc's. 



FIG. 2: (color online), (a) Differential conductance as a func- 
tion of magnetic field for a Pt-Ir/PrOs4Sbi2 point-contact 
junction at 90mK. (b) Field dependence of the zero-bias peak 
height and excess spectral defined in the text. 



sistive upper-critical field. Second, H"{T) and H'{T) 
appear to emerge from different low-temperature asymp- 
totes and gradually approach each other with increas- 
ing temperature. This non-parallel behavior between 
H"{T) and H'{T) resembles the H-T phase diagram 
determined from angular magneto-thermal conductivity 
in Ref.Q, in contrast to parallel phase boundaries de- 
termined from heat capacity in Ref.[37|. Finally, our 
H"{T) and H'{T) curves appear to approach zero at dif- 
ferent temperatures, wl.85K and «1.65K respectively, 
agreeing well with Td and Tc2 measured by x{T) on 
our samples. Here it is important to emphasize that 
our H-T diagram, which indicates two distinct super- 
conducting OPs with different symmetries and Tc's, is 
determined spectroscopically with a highly local probe, 
thus arguing strongly against sample inhomogeneity as 
the cause of multiple Tc's. It should also be noted that 
although our H-T diagram is qualitatively similar with 
the one reported in Ref.Q, there is quantitative dif- 
ference in the low-temperature asymptote between our 
H'{T) and their lower- field boundary H*{T). Interest- 
ingly, surface impedance measurements on PrOs4Sbi2 



have also indicated a lower-field phase boundary with 
a zero-temperature asymptote of wl.ST 38], which is 
quantitatively consistent with our H'{T). 

To probe the effects of Ru-doping on the two supercon- 
ducting OPs observed, we also performed PCARS mea- 
surements on Pr(Oso.98Ruo.o2)4Sbi2 crystals. These 2% 
Ru-doped crystals have Tc « 1.8K and Hc2 ~ 2.2T at 
200niK. Figure EJa) plots the temperature evolution of 
a typical dl/dV spectrum. At 90mK a ZBP is also ob- 
served, but its height is lower than in the undoped case. 
Broad hump structures clearly emerge at w±1.4mV and 
«±0.6mV as indicated by S3 and 64 respectively in the 
figure. With increasing temperature, the ZBP decreases 
in height while the humps become narrower, and the 
spectrum completely flattens out above Tc « 1.8K. The 
inset of FigHJa) shows the spectral evolution with mag- 
netic field at 90mK. While the ZBP collapses into the 
hump structures above wl.ST, a spectral hump is still 
visible up to Ri2.2T, unlike the undoped case which shows 
no pronounced humps (FigO^a)). Figure [D^b) plots the 
temperature dependences of the ZBP height, the excess 
spectral area, and the positions of S3 and S4, using the 
same data reduction as for Figllfb). Comparing Fig|4fb) 
and FiglTJb), it is clear that while both cases show excess 
spectral area up to the resistive T^ FigUfb) shows a more 
rapid ZBP collapse with temperature. Using the two-OP 
Andreev reflection scenario from above, these observa- 
tions indicate that Ru-doping strengthens the non-nodal 
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FIG. 4: (color online). (a) Differential conduc- 
tance spectrum as a function of temperature for a Pt- 
Ir/Pr(Oso.98Ruo.o2)4Sbi2 point-contact junction in zero mag- 
netic field. The inset shows the spectral evolution versus field 
at 90 mK. (b) Temperature dependence of the ZBP height 
and excess spectral area, defined in the same way as in Fig.l. 



OP while weakening the nodal OP, in agreement with 
both heat-capacity and penetration-depth measurements 
0, which showed the rapid emergence of a non-nodal 
OP with Ru-doping. It is interesting to note that, for 
our PCARS data with 2% Ru-doping, the higher of the 
two temperature scales governing the spectral evolution 
appears to correspond to two hump structures, denoted 
by ^3 and ^4. These hump structures suggest that the 
non-nodal OP itself may have multiple sub-components. 
Further measurements of crystals with higher Ru-doping 
levels are under way to elucidate this possibility. 

In summary, we have performed point-contact Andreev 
reflection spectroscopy on PrOs4Sbi2 down to 90mK and 
up to STesla, and observed microscopic evidence for mul- 
tiple superconducting OPs. Detailed analysis of the spec- 
tral evolution revealed a H-T phase diagram character- 
ized by two distinct OPs, one of which clearly having gap 
nodes and the other one being non-nodal and becoming 
more robust with just 2% Ru-doping. Our observations 
corroborate similar results from bulk measurements, and 
highlight the unconventional nature of the pairing state 



in this heavy fermion superconductor. 
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